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XmODOCTXOV 

As  a  result  of  the  above-ground  nuelaar  teat  baa  treaty,  DMA 
researchers  have  ss— lined  relevant  non-nuclear  structured  environments  as 
si— lators  of  the  lato-tl—  nuclear  situation,  studies  of  these  natur- 
ally-ooemrr lag  eaviroo—ats  are  relevant  because  the  convective  instabil¬ 
ity  pros— a—  that  operate  la  hlgh-altltude  nuclear  phenomena  also  occur 
naturally .  of  the  naturally-occurring  disturbed  environments,  equatorial 
spread  r  is  associated  with  the  highest  levels  of  structure— particularly 
at  the  peak  of  the  solar  cycle,  lore in,  ee  have  analyzed  the  effects  of 
such  structure  —  the  traas ionospheric  propagation  of  radio  signals  from 
the  DU  wideband  satellite. 

Kleetroa-deaslty  irregularities  in  the  equatorial  ionosphere  provide 

considerable  information  about  the  processes  responsible  for  their 

generation.  Of  particular  Interest  to  cos— inleations  systems  is  the  shape 

of  the  irregularity  po— r  spectrum  over  spatial  scales  between  about  o.l 

and  10  km.  These  irregularities  have  been  measured  in  situ  by  satellite 

(Al-S,  83-4)  and  rockets  (Plunex,  Condor)  as  well  as  by  remote  means.  The 

spectra  frost  in  situ  measure— nts  are  relatively  sparse,  however,  compared 

with  a  much  more  readily  available  measurement  of  the  structure:  complex 

signal  scintillation.  Phase  scintillation  measured  over  translonospherlc 

paths  at  very-high  through  microwave  frequencies  provides  a  nearly  direct 

measure  of  integrated  Irregularity  structure,  at  least  at  weak  disturbance 
1" 

levels. 

During  the  operation  of  the  DMA  Wideband  satellite,  hundreds  of 
passes  of  equatorial  phase  scintillation  were  obtained  at  Ancon,  Peru,  and 
Kwajalein,  Marshall  Islands.  The  phase  spectral  measurements  made  using 


* 

References  are  listed  in  order  of  appearance  at  the  end  of  this  report. 


these  data  provided  much  of  tba  basis  for  DVA  pbanoaanology  modelling  of 
Irregularity  spectral  shape.  In  order  to  refine  that  spectral  model,  a 
selected  portion  of  tbe  Wideband  equatorial  data  were  recently  reproces¬ 
sed,  using  tbe  techniques  developed  for  the  HILAT  experiment.  Tbe  purpose 
of  the  study  was  to  obtain  a  characterization  of  the  phase  spectrum  under 
weak,  and  strong  scattering  conditions.  Varying  conditions  result  in  spec¬ 
tral  forms  that  are  very  different.  In  weak  scatter,  the  phase  spectrum 
reflects  the  in  situ  structure  while  the  strong  scatter  spectral  shape  is 
significantly  altered  by  diffraction. 


The  first  portion  of  this  report  is  concerned  with  weak  scatter  data 
only.  DHA  has  recently  accepted  a  two-component  irregularity  model,  and 
SKI  has  obtained  representative  values  of  the  slopes  and  break  frequency 
for  that  model  from  the  phase  data.  Fresnel  nulls  occur  in  most  of  the 
weak  scatter  phase  spectra;  the  spectral  fitting  must  carefully  avoid  that 
superimposed  structure.  In  Section  2,  we  outline  the  approach  we  have 
used  in  the  analysis,  and  present  the  observed  phase  spectral  charac¬ 
teristics. 

In  our  strong  scatter  analysis  presented  in  Section  3,  we  have  pur¬ 
sued  two  aspects  of  the  effects  of  diffraction  on  phase.  The  first  is  the 
alteration  of  the  temporal  phase  record  during  strong  scattering;  from 
this,  an  estimate  of  phase  error  as  a  function  of  intensity  disturbance 
has  been  obtained.  The  second  aspect  is  the  overall  shape  of  the  phase 
error  as  a  function  of  disturbance  level;  from  this,  some  broad 
conclusions  about  strong  scatter  Fresnel  effects  on  phase  can  be  drawn. 


SECTION  2 


WEAK  SCATTER  ANALYSIS 


2.1  ANALYSIS  RATIONALE. 

A  Major  concern  of  those  using  phase  spectra  to  Imply  irregularity 
structure  Is  that  the  particular  manifestations  of  diffraction  on  the 
phase  spectrum. are  unknown.  By  diffraction,  we  mean  any  departure  from 
the  nominal  linear  dependence  of  signal  phase  on  carried  frequency,  until 
the  diffraction  effects  can  be  specifically  identified,  any  conclusions 
regarding  phase  spectral  shape  should  be  based  on  weak  scatter  data  only. 

A  convenient  range  of  disturbance  level  to  work  with  which  satisfies  this 
criterion,  is  an  intensity  scintillation  level  of  0.2  <  s,  <  0.5.  The 
lower  limit  is  necessary  to  avoid  bias  of  the  results  caused  by  noise 
contamination.  The  upper  limit  follows  from  analysis  (which  is  presented 
in  Section  3);  above  this  level  of  S^.  phase  errors  begin  to  rapidly 
increase  and  the  phase  spectral  shape  begins  to  be  distorted.  For  the 
passes  we  have  selected,  virtually  all  of  the  equatorial  Wideband  data 
that  fall  into  this  weak  scatter  disturbance  range  are  at  413  MHz  and  1239 
MHz. 


Limiting  the  weak  scatter  analysis  to  low  values  of  s„  is  not  with- 
out  its  complications.  In  particular,  Fresnel  null  structure  is  super¬ 
imposed  on  most  of  the  weak  scatter  phase  spectra  we  observe.  At  413  MHz, 
the  weak  scatter  often  arises  from  bottomslde  spread-F  irregularities 
which  are  contained  in  relatively  narrow  and  homogeneous  layers  below  the 
F- layer  peak.  For  Fresnel  structure  to  occur  in  the  phase  spectra,  the 
irregularity  layer  must  fit  the  thin  screen  definition.  Bottomside  spread 
F  approximates  a  thin  screen,  and  Fresnel  nulls  should  be  conspicuous.  At 
1239  MHz,  the  irregularities  producing  weak  scatter  scintillations  are 
much  stronger.  At  the  local  time  sampled  by  the  Wideband  satellite,  they 
are  often  the  decaying  remnants  of  plume  structures,  which  are  less  con- 


fined  in  altitude  and  lest  homogeneous  than  bottomside  structure.  Thus, 
we  expect  Fresnel  nulls  to  be  less  conspicuous  at  1239  KHz  than  at  413 
KHz.  and  this  is  the  case. 

Figure  1  shows  an  example  sequence  of  simultaneous  413  KHz  and  1239 
KHz  intensity  and  phase  spectra  showing  Fresnel  effects.  (Note  that  these 
1239  KHz  spectra  are  not  Included  in  our  weak,  scatter  results  because 
their  levels  were  less  than  0.2.)  The  arrows  indicate  the  expected 
frequency  of  occurrence  of  the  first  Fresnel  intensity  peak  and  phase 
null,  for  the  given  pass  geometry,  an  assumed  irregularity  anisotropy  of 
50:1  and  a  layer  height  of  350  km.  The  nulls  are  conspicuous,  even  in 
these  spectra  that  have  been  smoothed  for  display,  and  persist  even  up  to 

s.  levels  of  ~  0.7  at  413  KHz. 
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Unfortunately,  the  Fresnel  nulls  are  so  prevalent  in  the  413  mh t-  data 
that  they  complicate  the  extraction  of  spectral  slope  and  break  frequency 
information.  This  became  clear  when  a  multicomponent  fit  algorithm,  used 
early  in  the  processing,  gave  break  frequency  distributions  that  were 
dependent  upon  carrier  frequency.  This  tendency  contradicts  weak  scatter 
theory  assumption  that  a  break  in  the  irregularity  spectrum  should  appear 
in  the  phase  spectrum  at  the  same  temporal  frequency,  independent  of  the 
carrier.  Subsequent  numerical  simulations  using  the  multicomponent  fit 
algorithms  verified  that  the  carrier-dependent  break  distributions  were, 
in  fact,  due  to  the  Fresnel  nulls. 

To  avoid  such  Fresnel  null  biases,  we  chose  to  obtain  two  linear  fits 
over  temporal  frequency  ranges  near  the  low  and  high  frequency  ends  of  the 
phase  spectrum.  In  selecting  fit  ranges,  temporal  frequencies  affected  by 
the  primary  Fresnel  null  at  413  KHz  and  1239  KHz  were  avoided.  For  an  F- 
layer  altitude  of  350  km  and  an  irregularity  anisotropy  of  50:1,  the  413 
KHz  primary  null  is  centered  somewhere  between  0.7  Hz  and  1.3  Hz,  depen¬ 
ding  on  propagation  geometry.  For  1239  KHz,  the  range  is  between  1.2  Hz 
and  2.2  Hz.  Thus,  the  413  MHz  nulls  limit  the  fit  range  for  a  low  fre¬ 
quency  segment,  while  those  at  1239  KHz  limit  the  useable  high  frequency 
fit  range.  Taking  into  account  the  width  of  the  null,  convenient  fit 
ranges  are  0.10  Hz  to  0.56  Hz  (-1.00  to  -0.25  log  Hz,  approximately  4500  m 
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to  800  ■  spatial  scale),  and  1.8  Hz  to  10.0  Hz  (0.25  to  1.00  Hz,  approxi¬ 
mately  250  a  to  45  a  spatial  seals). 


Ths  outs Ids  raogs  Halts  of  0.1  and  10. 0  Hz  wars  chosen  prlaarily  to 

2 

maintain  historical  continuity  with  sarlisr  Wideband  analyses.  In  addi¬ 
tion.  however,  they  avoid  any  possible  low  frequency  spectral  contaaina- 
tion  caused  by  windowing,  and  high  frequency  contamination  caused  by 
systaa  noise.  Furtheraore,  if  the  power  spectral  density  at  10  Hz  was 
less  than  -50  dB  rad2/Hz  (a  conservative  estiaate  of  the  Wideband  system 
noise  floor),  the  high  frequency  fit  limit  was  chosen  to  end  at  that 
frequency. 

we  have  also  been  careful  to  avoid  contamination  of  the  results  by 
interference  or  by  other  spurious  effects.  This  was  done  by  applying 
relatively  stringent  requirements  to  each  set  of  simultaneous  137  MHz.  413 
HHz  and  1239  KHz  spectra.  Any  spectrua  that  was  accepted  had  an  inte¬ 
grated  energy,  which  was  consistent  with  that  at  the  other  two  frequen¬ 
cies.  More  specifically,  if  the  ratio  between  the  ras  phase  deviations  at 
two  adjacent  carriers  exceeded  30  percent  of  the  linear  ratio,  the  spectra 
were  rejected  frea  the  data  base. 

As  for  the  spectral  analysis,  per  se,  we  have  used  (up  through  the 
spectral  slope  fitting)  the  current  HXLAT  in-field  processing.  This 
Includes  edge-aatching  and  cosine-windowing  of  30-s  spans  of  data  at  125 
Hz  before  spectral  analysis.  Prior  to  obtaining  the  slopes,  the  phase 
power  spectrua  is  log-dec laa ted  to  50. 

For  the  reanalysls  of  the  Wideband  equatorial  data,  we  chose  24 
orbits  from  Ancon  and  19  from  Kwajaleln.  All  were  selected  because  they 
show  consistent  scintillation  throughout  aost  of  each  pass,  i.e.,  at  a 
variety  of  propagation  geometries.  The  Ancon  data  are  from  the  period 
February  through  April  1977;  the  Kwajaleln  data  are  froa  July  through 
September  1977.  These  are  the  peat  scintillation  seasons  at  each  longi¬ 
tude3  in  the  year  of  maximum  solar  activity  sailed  by  Wideband.  The 
spectral  parameter  statistics  froa  the  two  stations  are  virtually  iden¬ 
tical;  what  we  discuss  in  the  paragraphs  below  is  a  combined  result. 
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2.2  POWER-LAW  AHD  TVO-COMPONENT  POPULATIONS. 

Our  primary  interest  is  with  two-component  phase  spectra  for  which 
the  low  frequency  slope.  P  ,  is  shallower  than  the  high  frequency  slope, 

A 

Pb«  This  spectral  form  does,  in  fact,  comprise  the  majority  of  our  obser¬ 
ved  spectra.  There  is  also  a  population  of  spectra  for  which  the  two 
slopes,  P_  and  P.  ,  are  very  close  to  one  another;  as  a  practical  matter, 

a  D 

these  are  single  power-law  in  form.  The  third  category  of  spectra  are 
those  for  which  P_  is  significantly  steeper  than  P.  .  Although  approxi- 
mately  30  percent  of  the  spectra  are  of  this  form,  we  have  discarded  them 
from  the  statistics.  Our  justification  is  that  this  spectral  form  is 
almost. never  observed  in  situ  at  the  equator.  Host  likely,  the  phase 
spectra  of  this  type  are  a  result  of  integration  through  multiple  layers 
having  different  characteristics. 

The  working  data  set  consists,  then,  of  the  two-component  spectra 
with  P_  <  P.  ,  and  the  P  -  P.  single  power- law  cases.  Spectra  are 

AD  AD 

identified  as  single  power-law  if  P  is  within  20  percent  of  P.  .  This 

A  D 

criterion  is  somewhat  arbitrary;  the  specific  percentage  chosen  does 
affect  the  slngle/dual  power-law  population  mix,  but  does  not  significant¬ 
ly  alter  the  slope  or  spectral  break  distributions.  Using  this  criter¬ 
ion,  the  percentage  of  single  power-law  spectra  is  37  at  413  MHz,  and  28 
at  1239  MHz.  These  percentages  are  not  a  strong  function  of  carrier  fre¬ 
quency,  suggesting  that  the  same  irregularity  process  or  mix  of  processes 
is  operating  over  a  wide  range  of  turbulence  strength.  For  example,  the 
two-component  spectra  do  not  appear  to  be  limited  to  the  strong  turbulence 
that  produces  1239  MHz  scintillation. 

2.3  SPECTRAL  SLOPES. 

The  distributions  of  the  low-  and  high-frequency  slopes,  P  and  P.  , 

A  D 

are  shown  in  Figure  2  for  the  combined  Ancon  and  Kwajaleln  spectra  (ex¬ 
cluding  those  for  which  P  is  much  steeper  than  P.  ,  as  noted  above).  The 

A  D 

data  for  413  MHz  and  1239  MHz  are  shown.  The  weak  scatter  population  at 
137  MHz  is  too  small  to  provide  useful  distributions,  although  the  137  MHz 
mean  values  of  P_  and  P.  are  similar  to  those  at  the  higher  frequencies. 

A  D 


The  sample  populations  and  tha  naan  valua  for  each  distribution  are  noted. 
The  shapes  of  tha  distributions  for  P_  and  P.  are  similar  for  the  two 
carrier  frequencies,  and  their  naans  are  nearly  the  same,  with  P  ~  2.3 

cl 

and  P.  *•  3.9. 

D 

The  agreement  between  the  two  carriers  is  an  important  result  because 
it  shows  that  (in  weak  scatter,  at  least)  there  is  no  appreciable  shallow¬ 
ing  of  the  spectral  slope  of  either  component  with  disturbance  level.  As 
mentioned  previously,  the  irregularity  structure  that  produces  weak  scat¬ 
ter  scintillation  at  1239  KHz  is  several  times  stronger  than  that  causing 

2.3 

the  same  signal  disturbance  at  413  KHz.  Previous  studies  and  implica¬ 
tions  strong  scatter  temporal  statistics5  showed  a  shallowing  of  slope 
with  increasing  disturbance  level.  There  should  be  explanations  for  this 
apparent  discrepancy.  In  the  case  of  Reference  2,  the  difference  is 
caused  by  diffraction  effects,  as  the  data  in  Section  3  will  show. 

2.4  SPECTRAL  BREAK  SCALES. 

In  addition  to  the  slopes,  the  spatial  scale  at  which  the  two- 
component  spectrum  breaks  is  an  Important  factor  in  characterizing  the 
Irregularity  generation  process.  For  phenomenological  purposes,  the 
spatial  size  at  which  this  break  occurs  is  termed  the  freezing  scale.  In 
the  propagation  case,  a  break  in  the  irregularity  spectrum  appears  in  the 
phase  spectrum  at  a  temporal  frequency  dictated  by  the  line-of-sight  scan 
velocity.  This,  in  turn,  depends  upon  the  Irregularity  altitude, 
anisotropy,  and  motion.  The  temporal  frequency  at  which  the  translated 
break  occurs  is  independent  of  carrier  frequency,  as  previously  noted. 

This  lends  credence  to  the  Interpretation  that  its  origin  is  in  situ. 

The  orbital  Inclination  of  the  Wideband  satellite  is  such  that  for 
the  premldnlght  pass  data  used  here,  the  propagation  path  scans  nearly 
along  the  magnetic  meridian.  For  F-region  irregularities  that  are  highly 
extended  along  the  magnetic  field  (50:1  anisotropy),  the  effective  scan 
velocity  caused  by  the  satellite  motion  has  an  average  value  of  ~  0.45 
km/s.  At  the  local  time  and  latitude  of  the  observations,  F-region 
polarization  and  tidal  electric  fields  will  typically  drive  the  irregu- 
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larities  eastward  at  a  velocity  of  ~  0.25  ka/s.  The  total  effective  scan 
velocity  will  then  be  in  the  0.5  ka/s  to  0.7  ka/s  range. 


For  each  spectrua.  the  break  frequency  is  taken  to  be  at  the  inter¬ 
section  of  the  P.  and  P.  linear  fit  segaents.  The  distribution  of  the 
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break  frequencies  is  shown  in  Figure  3  for  413  KHz  and  1239  KHz.  The 
distributions  for  both  carriers  are  broad,  and  are  slallar  in  shape  and 
mean  value.  Much  of  the  broadening  of  the  distribution  is  because  the 
scan  velocity  of  the  propagation  path  through  the  irregularities  depends 
upon  the  propagation  geometry,  a  break  at  a  fixed  spatial  scale  will 
therefore  appear  at  different  temporal  frequencies  throughout  the  course 
of  a  pass.  If  there  were  a. means  of  accurately  applying  correction  for 
the  scan  velocity  so  that  the  data  could  be  sorted  in  spatial  scale,  the 


distributions  would  be  much  narrower. 
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Figure  3.  Distributions  of  phase  spectral  break  frequencies. 
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Choosing  an  overall  average  break,  frequency  of  1.07  Hz.  and  a  repre¬ 
sentative  scan  velocity  of  0.60  ka/s,  our  estimate  of  the  irregularity 
break  scale  is  near  0.65  ka.  This  is  in  good  agreement  with  the  irregu¬ 
larity  spectrua  break  scale  observed  by  satellite  and  rockets — typically 
reported  as  0.75  ka.  The  difference  is  probably  not  significant,  because 
the  phase-implied  value  is  highly  dependent  upon  our  scan  velocity  assump¬ 
tions  about  layer  height  and  irregularity  anisotropy. 

2.5  CO HP AS IS OH  WITH  PREVIOUS  STUDIES. 

It  is  worthwhile  to  compare  the  original  Wideband  slope  results1  with 
those  obtained  here.  For  the  earlier  Wideband  analysis,  the  slopes  were 
extracted  between  0.1  Hz  and  10.0  Hz  (as  has  been  done  here),  but  it  was 
presumed  that  the  spectra  were  of  single  power-law  fora.  There  were 
additional  differences,  however,  in  the  pre-fit  processing,  which  could 
have  bad  an  effect  on  the  spectral  index  statistics.  These  were:  (1) 
prior  to  spectral  analysis,  the  temporal  records  were  filter  detrended  but 
not  windowed,  and  (2)  no  requirements  were  placed  on  rms  phase  deviation 
ratios,  nor  were  precautions  taken  to  avoid  noise  contamination  of  the 
fits  at  low  disturbance  levels. 

To  determine  whether  the  current  pre-fit  reprocessing  (windowing  and 
stringent  sorting)  has  significantly  affected  the  slope  statistics,  we 
have  made  single  power-law  fits  (0.1  Hz  to  10.0  HZ)  for  our  reprocessed 
passes.  The  results  are  shown  in  Figure  4,  adjacent  to  the  replotted 
results.  We  can  see  that  the  current  analysis  slopes  are  systematically 
steeper  than  those  from  the  original  analysis.  For  levels  less  than 
about  0.3.  this  is  not  unexpected  because  the  Wideband  spectral  fit  range 
included  the  noise  floor  [Item  (2),  first  paragraph  of  this  section]. 
Reference  1  noted  this  fact,  and  in  our  reanalysis,  we  have  specifically 
avoided  noise  contamination. 

For  >  0.3,  where  noise  is  not  a  problem  in  the  spectral  fits,  the 
revised  slopes  are  also  steeper  than  the  original  by  about  5  to  10  per¬ 
cent.  This  is  most  likely  an  effect  of  sidelobe  contamination  of  the  high 
frequency  portion  of  the  original  spectra.  This  is  a  common  problem  with 
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spectral  analysis  of  power -law  processes;  we  now  eliminate  the  sidelobes 
by  windowing  the  temporal  records. 

The  revised  value  of  slope  near  3.0  is  undoubtedly  a  more  accurate 
representation  of  the  underlying  irregularity  spectrum  than  were  the 
original  Wideband  values  of  -  2.5.  The  difference  does  point  out.  how¬ 
ever.  that  there  was  a  deficiency  in  the  earlier  Wideband  processing,  the 
difference  is  of  little  consequence.  In  terms  of  understanding  irregular¬ 
ity  processes,  a  much  more  important  shortcoming  of  the  original  Wideband 
analysis  was  the  presumption  that  the  phase  spectra  were  purely  power-law. 


SECTION  3 


STRONG  SCATTER  AH  ALTS  IS 

The  m)i  scatter  pitas*  spectral  fora  remains  consistent  up  through  a 
disturbance  level  for  which  the  index  is  approxiaately  0.5.  Above  this 
level,  however,  the  shape  of  the  phase  spectrum  is  significantly  and  pro¬ 
gressively  altered  by  diffraction  processes.  As  an  example.  Figure  5 
shows  this  change  of  spectral  shape  in  terms  of  the  mean  value  of  P  and 
Pfe  indices  at  137  KHz  and  413  KHz  during  moderate  to  strong  scatter.  The 
slopes  of  both  the  low  and  high  frequency  portions  of  the  spectra  depart 
from  the  weak,  scatter  values  of  P.  -  2.3  and  P.  -  3.9  for  S,  >  0.5,  and 
systematically  become  nor*  shallow  at  stronger  scatter  levels. 

we  will  consider  two  aspects  of  the  behavior  of  phase  scintillation 
under  strong  scatter  conditions.  The  first  is  a  look  at  the  details  of 
phase  temporal  behavior  in  strong  scatter.  This  provides  estimates  of 
diffraction-induced  phase  errors  and  how  they  arise,  we  then  give  an 
overall  view  of  the  phase  spectral  shape  as  a  function  of  scattering 
level,  from  which  some  broad  conclusions  regarding  Fresnel  dependence  can 
be  made. 

3.1  TEMPORAL  CHARACTER  OF  PHASE  IH  STRONG  SCATTER. 

The  Wideband  experiment  provided  a  nearly  ideal  means  to  study  the 
onset  and  development  of  diffraction  effects  on  signal  phase.  This  was 
the  Wideband  UHF  comb,  seven  carriers  transmitted  by  the  spacecraft  beacon 
at  a  spacing  of  approxiaately  11  MHz,  between  378  MHz  and  435  MHz.  We 
have  used  these  data  to  estimate  the  nature  and  degree  of  alteration  of 
the  phase  temporal  record  by  diffraction.  The  reasoning  is  that  two  UHF 
signals  separated  by  only  11  MHz  should  undergo  nearly  identical  propaga¬ 
tion  effects,  while  those  separated  at  multiples  of  11  MHz  up  through  66 
MHz  should  exhibit  increasingly  different  diffraction  effects.  It  is  thus 


possible  to  observe  diffraction-induced  detail  in  the  signal  structure, 
which  would  not  be  possible  In  a  comparison  between  more  widely  spaced 
carriers. 
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For  the  analysis,  we  chose  a  high-elevation  Wideband  pass  from  Ancon, 
for  which  the  UHF  scintillation  ranged  from  weak,  to  saturated.  The  first- 
order  statistics  (S4  indices  and  ras  phase  deviations)  for  the  137  MHz  and 
413  MHz  wideband  carriers  are  shown  in  Figure  6.  As  can  be  seen,  the 
scintillation  change  is  relatively  systematic  over  the  duration  of  the 
pass  and  shows  no  particular  dependence  on  propagation  geometry. 

Figure  7  shows  the  signal  intensity  and  phase  records  for  one  30-s 
portion  of  the  pass,  comparing  the  413  MHz  and  378  MHz  carriers.  The 
phase  data  have  been  edge-matched  and  windowed,  i.e.,  as  prepared  for 
routine  spectral  analysis.  Also  shown  is  the  difference  between  the  phase 
records  at  the  two  carriers,  following  the  slight  but  necessary  linear 
scaling  for  measurement  frequency. 

For  this  particular  span  in  the  figure,  the  scintillation  is  moderate 
(S4  -  0.6),  and  even  though  the  carriers  are  separated  by  33  MHz.  the 
phase  records  are  virtually  indistinguishable.  This  is  confirmed  by  the 
phase  difference  plot,  which  is  near  zero,  although  a  series  of  small  and 
rapid  excursions  can  be  seen.  These  are  caused  by  cusp-like  changes  in 
phase,  which  occur  in  the  signal  minima  of  intensity  fades.  In  this  case, 
where  the  fades  are  only  a  few  decibels,  the  cusps  have  deviations  that 

are  only  a  few  degrees  from  the  nominal  phase  record. 

• 

Figure  8  shows  a  second  30-s  span,  where  S4  -  0.75.  Diffraction  has 
produced  quite  different  intensity  records  at  the  two  carriers,  and  the 
corresponding  phase  cusps  are  large  enough  to  significantly  alter  the 
phase  record.  The  two  records  are  similar  only  in  their  large-scale 
variation;  their  difference  is  step-like,  due  to  the  large  phase  changes 
that  occur  during  fades  at  one  or  the  other  carrier. 

Finally,  Figure  9  is  a  third  30-s  span  for  near-saturated  scintilla¬ 
tion  (S4  -  0.9).  The  near-continuous  fading  of  both  carriers  produces 
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phase  records  that  are  not  even  similar,  thus  the  phase  difference  is  more 
random  and  has  lost  the  step-like  character  of  the  previous  figure. 

From  these  examples,  it  can  be  seen  that  phase  cusps  during  fades  are 
the  primary  source  of  diffraction  contamination  in  the  phase  spectra  for 
moderate  through  strong  scatter.  A  key  point  is  that  the  cusps  in  each 
phase  record  depend  upon  the  details  of  the  diffraction-generated  signal 
intensity  structure,  and  those  details  change  very  rapidly  with  measure¬ 
ment  frequency.  As  a  result,  even  for  closely  spaced  carriers,  minor 
differences  in  the  characteristics  of  the  two  fading  signals  can  result  in 
reconstruction  of  very  different  temporal  phase  records,  we  point  out 
that  intensity  correlation  between  carriers,  per  the  conventional  defini¬ 
tion5  reflects  very  little  of  the  Intensity  structure  detail  that  controls 
the  phase  cusps;  for  example,  the  intensity  records  in  Figure  9  are  cor¬ 
related  above  80  percent. 

It  is  also  useful  to  note  that  the  phase  cusp  behavior  during  fades 
is  not  dependent  (to  first  order)  upon  data  sampling  rate.  The  phase  data 
in  Figures  7  through  9  were  constructed  from  quadrature  components  using  a 
sampling  rate  of  125  Hz.  We  have  repeated  the  process  using  the  full 
Wideband  sampling  rate  of  500  Hz,  and  have  obtained  records  nearly  identi¬ 
cal  to  those  from  the  lower  reconstruction  rate.  To  see  no  change  despite 
a  four-fold  increase  in  temporal  resolution  shows  that  the  phase  cusps  are 
adequately  sampled,  and  are  not  caused  by  a  failure  to  keep  track  of 

rapidly  changing  phase  during  a  fade. 
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In  Figure  10,  the  standard  deviation  of  the  diffraction-induced  phase 
differences  are  plotted  as  a  function  of  s^.  Each  point  corresponds  to 
30-s  spans  such  as  those  shown  in  Figures  7  through  9.  Three  carrier 
separations  are  shown,  66  MHz,  33  MHz,  and  11  MHz;  the  66  MHz  spacing  data 
are  from  the  outside  carrier  pair  of  the  UHF  comb  (435  MHz  and  378  MHz); 
the  33  KHz  and  11  MHz  spacing  data  are  from  various  carrier  pair  combina¬ 
tions  within  the  comb.  A  least-squares  quadratic  fit  has  been  made  to 
each  block  of  data  to  illustrate  that  the  diffraction-induced  phase  dif¬ 
ference  dependence  on  index  is  nearly  the  same  for  all  three  carrier 
separations.  We  note  that-  in  weak  scatter,  the  differences  are  small  (~ 


20  deg),  even  for  the  66  MHz  spacing;  this  is  most  likely  not  due  to  phase 
cusps.  Zt  is  apparently  a  result  of  differential  Fresnel  null  structure. 
In  addition,  there  is  dispersion  in  the  spacecraft  antenna,  which  produces 
the  weak  large-scale  trend  that  can  be  seen  in  Figure  7. 

The  dependence  of  the  diffraction-induced  phase  errors  on  carrier 
separation  are  shown  more  clearly  in  Figure  10(d),  where  the  least-squares 
fits  are  compared.  As  the  level  of  disturbance  increases  from  weak 
scatter  to  S4  -  0.6,  the  separation  of  the  curves  remains  nearly  constant. 
The  curves  then  begin  to  converge  until  at  S4  ~  0.8,  the  phase  differences 
become  independent  of  carrier  separation  at  a  value  near  8  rad. 

The  important  point  to  be  made  from  Figure  10  is  that  diffraction  can 
create  significant  differences  in  phase  structure  even  between  closely 
spaced  measurement  frequencies.  Thus,  even  at  moderate  scatter  levels  (S4 
~  0.6),  there  will  be  differences  in  phase  spectral  shapes  measured  at  UHF 
carriers  separated  by  only  a  few  tens  of  megahertz.  Again,  this  is 
because  the  phase  structure  depends  almost  entirely  on  the  details  of  the 
diffraction-developed  intensity  structure. 

It  is  possible  to  use  the  result  in  Figure  10  to  estimate  the  dif¬ 
fraction  contamination  that  would  be  encountered  in  the  case  of  a  single 
carrier  measurement.  First,  we  note  that  the  rms  phase  differences  in  the 
figure  include  contamination  from  phase  cusps  at  both  carriers.  As  the 
carrier  separation  increases,  so  does  the  dissimilarity  of  the  intensity 
fade  patterns;  thus,  the  total  number  of  phase  cusps  from  the  two  carriers 
increases.  In  this  sense,  the  carrier  separation  dependence  of  the  phase 
difference  in  Figure  10  is  a  product  of  (and  is  applicable  only  to)  the 
closely  spaced  carrier  measurement.  In  the  case  of  the  single  carrier 
measurement,  diffraction  structure  at  any  frequency  would  be  no  different 
from  that  at  any  one  of  the  comb  carriers  for  the  same  disturbance  level: 
Yet,  single  carrier  phase  would  be  compared  to  a  coherent  reference  phase 
from  a  carrier  at  least  1000  MHz  away,  and  presumably  undisturbed.  Thus, 
the  single  carrier  diffraction  effect  in  phase  should  be  comparable  to  or 


less  than  those  from  our  most  closely  spaced  comb  carriers,  i.e.,  a  few 

degrees  at  S.  ~  0.5  and  a  few  radians  at  S.  -  0.9. 
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This  result  is  also  in  satisfactory  agreement  with  the  average  spec¬ 
tral  shape  data  shown  in  Figure  5.  At  -  0.5,  where  the  integrated 
phase  error  due  to  diffraction  is  only  a  few  degrees,  there  can  be  little 
spectral  shape  distortion,  as  is  the  case.  On  the  other  hand,  when  the 
scintillation  is  saturated,  we  expect,  and  see,  significant  distortion  in 
both  the  high  and  low  frequency  portions  of  the  phase  spectrum  as  is  con¬ 
sistent  with  an  rms  phase  differences  of  a  few  radians. 

3.2  PHASE  AND  INTEHSITY  SPECTRAL  SHAPES  IN  STRONG  SCATTER. 

Knowing  something  of  the  manifestations  of  the  diffraction  process  in 
the  phase  temporal  records,  we  now  turn  to  a  brief  review  of  how  that 
structure  is  distributed  in  the  phase  spectrum,  i.e.,  the  spectral  shape 
in  strong  scatter.  Figure  5  provides  a  glimpse  of  the  spectral  shape 
dependence  on  disturbance  levels  In  terms  of  P_  and  P.  indices,  which 
shallow  at  both  low  and  high  frequencies.  In  strong  scatter,  of  course, 
characterization  of  the  spectrum  by  and  is  probably  inappropriate, 
because  the  spectral  form  may  be  altered  or  disguised  by  diffraction 
effects . 

Strong  scatter  theories  abound  in  the  literature,  and  there  are 
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various  approaches  to  the  problem.  *  '  *  By  and  large,  these  theories 

and  numerical  simulations  address  only  the  signal  intensity  behavior,  and 
most  are  limited  to  gausslan  or  power-law  source  spectral  shapes.  There 
is  generally  good  agreement  between  observation  and  predicted  strong  scat¬ 
ter  intensity  spectra.  The  shape  of  the  intensity  spectrum  in  strong 
scatter  is  dominated  by  the  Fresnel  scale  for  the  particular  propagation 
geometry  and  operating  frequency,  just  as  it  is  in  weak  scatter.  As  the 
perturbation  level  increases  from  weak  through  strong  scatter,  the  narrow 
peak  of  the  weak  scatter  intensity  spectrum  broadens  toward  both  high  and 
low  frequencies.  This  bulge  has  the  effect  of  shallowing  the  low  frequen¬ 
cy  portion  of  the  spectrum,  but  steepening  the  high  frequency  end,  because 
the  Intensity  structure  at  short  spatial  scales  remains  constant.  Further- 


more,  between  the  Fresnel  scale  and  the  high  frequency  limit  of  the  inten¬ 
sity  structure,  the  spectrum  is  close  to  gaussian  in  shape,  having  con¬ 
tinuous  curvature.  This  means  that  because  of  the  Fresnel  scale  depen¬ 
dence  on  carrier  frequency,  spectral  slopes  extracted  from  the  same  high 
frequency  portion  of  Intensity  spectra  will  systematically  depend  upon 
that  frequency. 

Very  little  work  has  been  done  with  strong  scatter  phase  behavior, 
either  observationally  because  of  measurement  difficulty,  or  theoreti¬ 
cally  because  of  inadequacies  in  the  propagation  models.  The  conventional 
opinion  is  that  the  phase  spectrum  should  become  more  like  the  intensity 
spectrum  as  the  scintillation  level  approaches  saturation.10  Presumably 
then,  we  should  see  a  transition  from  the  weak  scatter  two-component  form 
to  one  with  a  broadening  of  energy  centered  on  the  Fresnel  scale. 

The  reprocessed  Wideband  data  affords  an  excellent  opportunity  to 
investigate  strong  scatter  effects  on  the  phase  spectrum.  However, 
because  of  the  uncertainty  about  distortion  of  the  spectral  shape  by 
diffraction,  we  have  chosen  not  to  extract  spectral  indices  by  curve 
fitting.  Instead,  the  spectra  have  been  sorted  into  ranges  and 
averaged.  This  approach  is  not  without  its  shortcomings;  in  particular, 
the  average  spectrum  tends  to  be  dominated  by  the  most  strongly  disturbed 
spectra  in  each  range.  The  essence  of  the  spectral  shape  is  not  lost, 
however,  as  is  shown  in  Figure  11,  where  the  average  phase  and  intensity 
spectra  are  shown  for  the  weak  scatter  (S^  range  0.3  to  0.5)  at  413  MHz. 
Fits  to  this  averaged  phase  spectrum  over  the  selected  frequency  ranges 
for  P  and  P  give  indices  quite  close  to  those  in  the  slope  distributions 
of  Figure  6;  P.  is  2.3,  as  was  found  from  individual  fits,  while  Pw  is  3.7 
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compares  well  with  the  distribution  value  of  3.9.  Thus,  we  feel  that 
spectral  forms  obtained  from  averaging  are  accurate  enough  to  be  a  useful 
means  to  characterize  the  complex  signal  behavior  in  strong  scatter. 

In  passing,  two  additional  points  can  be  made  from  the  weak  scatter 
spectrum  in  Figure  11.  First,  even  after  averaging  over  a  wide  variety  of 
propagation  geometries,  the  spectrum  still  clearly  shows  the  Fresnel  null 
near  1  Hz;  this  is  an  indication  of  how  prevalent  that  propagation  feature 
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Figure  11.  Average  weak  scatter  spectral  forms. 

Is  In  weak  scatter  spectra.  A  second  point  Is  that  at  first  glance,  the 

spectrum  appears  to  be  of  single  power-law  form — even  though  we  know  it  is 

a  component  spectrum  having  very  different  values  of  P  and  P.  .  This 

a  b 

emphasizes  the  need  for  consistent  and  careful  numerical  (not  visual) 
curve  fitting  In  comparative  studies. 

Figure  12  is  a  comparison  sequence  of  weak  through  saturated  scintil¬ 
lation  Intensity  and  phase  spectra  at  413  KHz.  The  S4  ranges  and  sample 
populations  for  these  average  spectra  are  as  indicated.  For  reference  on 
the  phase  spectra,  the  overall  weak  scatter  P  and  P.  regimes  are  shown. 

a  D 

The  Intensity  spectra  in  the  figure  generally  exhibit  the  behavior  predic¬ 
ted  by  strong  scatter  theory.  Starting  from  the  weak  scatter  form,  where 
the  intensity  spectrum  peaks  at  an  average  Fresnel  scale,  the  spectra 
broaden  as  the  perturbation  strength  increases.  As  noted  above,  the 
effect  is  a  bulging  of  the  spectrum  from  a  point  near  the  Fresnel  scale  to 
both  lower  and  higher  frequencies. 
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Figure  12.  Average  strong  scatter  spectral  forms  as  a  function  of  disturbance  level. 

The  corresponding  changes  In  the  phase  spectra  between  weak,  and 
strong  scatter  are  less  conspicuous,  except  for  a  general  upward  shift  as 
the  phase  variance  Increases.  However,  the  shallowing  of  the  overall 
spectrum  Indicated  by  Figure  5,  can  be  seen.  Figure  12  also  shows  that 
the  Increase  In  phase  power  spectral  density  Is  most  rapid  at  higher 
temporal  frequencies. 

a 

In  order  to  look  more  closely  at  the  change  in  spectral  shape,  and  to 
see.  In  particular,  if  there  are  Fresnel  dependences  in  phase,  differen¬ 
tial  power  spectral  densities  have  been  calculated  for  different  S4 
values.  These  are  shown  in  Figure  13  for  137  MHz  and  413  MHz.  The 
quantities  shown  are  the  differences  between  the  phase  power  spectral 
density  at  the  S4  range  indicated  and  those  In  the  average  spectrum  for  an 
S4  range  of  0.4  to  0.6.  The  Intensity  plots  show  how  consistently  the 
Fresnel  scale  dominates  the  Intensity  spectral  shape.  At  the  Fresnel 
scale,  the  power  spectral  density  remains  nearly  constant,  even  Into 
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Figure  13.  Differential  power  spectral  density  as  a  function  of  disturbance  level, 
referenced  to  average  power  spectral  density  for  0.4  <  S4  <  0.6. 
Relative  Fresnel  scales  are  indicated. 


saturation.  Lines  indicate  the  expected  relative  locations  of  the  137  MHz 
and  413  MHz  Fresnel  scales,  and  they  line  up  precisely  with  nulls  in  the 
differential  spectral  densities.  The  figure  also  shows  that  there  is 
little  change  in  the  frequency  at  which  the  Intensity  spectrum  changes 
slope,  despite  the  large  change  in  spectral  energy. 

The  differential  power  spectral  densities  also  show  that  the  phase 
spectra  in  strong  scatter  have  a  less  obvious  but  recognizable  dependence 
on  Fresnel  scale.  If  we  first  consider  the  413  MHz  data  [Figure  13(b)] 
there  are  two  regimes  in  differential  spectral  density,  separated  by  a 
peak-null  pair  near  1  Hz.  The  peak,  of  this  feature  is  aligned  with  the 
Fresnel  scale,  and  is  just  the  primary  phase  null  that  was  our  concern  in 
the  weak  scatter  analysis.  At  spatial  scales  larger  than  the  Fresnel 
scale,  the  spectral  energy  gained  with  increasing  disturbance  is  more 
uniform  and  less  rapid  than  it  is  at  shorter  spatial  scales. 

At  137  MHz,  the  Fresnel  null  effect  in  the  differential  spectral 
density  is  also  clear.  Hhat  is  less  clear,  relative  to  the  413  MHz  data, 
is  that  there  is  a  distinct  demarcation  between  the  rate  of  energy  gain  at 
spatial  scales  above  and  below  the  Fresnel  scale.  If  the  scale  size 
dependence,  which  is  so  clear  in  the  413  MHz  data,  is  purely  a  propagation 
effect,  it  should  be  observed  at  137  KHZ. 

Ve  speculate  that  the  difference  in  scale  size  dependence  between  the 
differential  power  spectral  densities  at  137  MHz  and  413  MHz  is  a  result 
of  Fresnel  effects  acting  on  a  two-component  irregularity  spectrum.  The 
nominal  freezing  scale  in  the  two-component  spectrum  is  near  750  m,  which 
is  very  near  the  413  KHz  Fresnel  scale  for  the  average  Wideband  propaga¬ 
tion  geometry.  At  137  KHz,  the  corresponding  Fresnel  scale  is  near  1300 
m.  In  the  case  where  the  Fresnel  and  freezing  scales  are  aligned,  the 
strong  scatter  spectrum  could  reflect  and  even  accentuate  the  underlying 

9 

two-component  nature  of  the  source  spectrum.  Simulations  have  shown  that 
the  strong  scatter  intensity  spectrum  depends  critically  on  the  slope  of 
the  source  spectrum  in  the  vicinity  of  the  Fresnel  scale;  this  may  be  even 
more  the  case  with  phase. 


A  more  detailed  analysis  of  the  dependencies  of  the  phase  spectra  is 
beyond  the  scope  of  this  report.  Independent  of  the  interpretation  of  the 
spectral  shape,  these  data  provide  a  precise  measure  of  the  strong  scatter 
effects  on  phase  at  frequencies  of  interest  to  communications.  Further¬ 
more,  they  present  guidelines  with  which  to  test  the  validity  of  future 
strong  scatter  theory  and  simulations  of  phase  behavior. 
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SECTIOI  4 
CONCLUSIONS 

We  have  addressed  the  question  of  phase  spectral  shape  in  the  10  tan 
to  0.1  tan  spatial  scale  regies,  under  weak,  and  strong  scattering  condi¬ 
tions.  He  used  nmltif requency  complex  signal  data  from  equatorial  obser¬ 
vations  of  the  ONA  Wideband  satellite,  that  we  have  reprocessed  using 
current  HXLAT  analysis  techniques.  The  purpose  of  the  study  was  twofold: 
(1)  to  characterize  the  two-component  phase  spectral  form  in  weak  scatter, 
which  can  be  simply  related  to  in  situ  irregularity  structure,  and  (2)  to 
obtain  estimates  of  the  effects  of  diffraction  on  the  phase  spectrum 
during  strong  scatter. 

He  have  found  that  a  useful  weak  scatter  range  corresponds  to  inten¬ 
sity  scintillation  levels  up  through  S4  -  0.5.  Up  to  this  level,  diffrac¬ 
tion  errors  are  much  less  than  a  radian,  and  the  phase  spectral  shapes  are 

consistent  and  independent  of  carrier  frequency.  The  average  two-compo- 

—2  3 

nmnt  phase  spectrum  has  a  low  frequency  dependence  near  f  and  a  high 

—3  9 

frequency  dependence  near  f  *  ,  separated  by  a  break  in  the  spectrum  near 
650  m.  Two-component  and  power-law  spectra  comprise  the  bulk  of  the  weak 
scatter  data;  of  these,  the  two-component  form  is  more  prevalent.  Host  of 
the  weak  scatter  data  at  UHF  show  Fresnel  nulls,  and  it  is  essential  that 
these  be  avoided  in  order  to  obtain  accurate  measures  of  the  in  situ 
spectral  indices. 

From  our  strong  scatter  analysis,  we  have  estimated  the  effects  of 
diffraction  on  the  phase  signal  and  on  the  phase  spectrum.  Using  closely 
spaced  carriers  from  Wideband,  the  nature  of  diffraction  effects  on  the 
phase  signal  have  been  determined.  Diffraction  contamination  for  any 
signal  depends  most  critically  upon  the  minor  details  of  its  intensity 
fading  pattern.  Quantitative  estimates  of  diffraction  contamination  show 
that  even  between  closely  spaced  carriers,  differential  structure  will 
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occur  at  moderate  scintillation  levels.  As  for  the  overall  effects  of 


diffraction  on  the  phase  spectrum,  we  have  obtained  average  strong  scatter 


spectral  shapes  as  a  function  disturbance  level.  There  is  evidence  in 


these  data  of  a  carrier  dependence  of  phase  spectral  form.  We  attribute 
this  dependence  to  the  coincidence  of  Fresnel  and  freezing  scales  at  UHF, 
and  their  separation  at  VHF. 


Ultimately,  the  propagation,  results  presented  here  must  be  used  to 
make  predictions/extrapolations  to  the  nuclear  disturbed  environment. 

This  involves  some  key  assumptions  that  must  be  addressed  in  future  work. 
One  is  the  degree  to  which  the  phase  spectral  shape  directly  reflects  in 
situ  structure.  This  report,  for  the  first  time,  has  addressed  questions 


of  diffraction  contamination  as  a  function  of  s,  value.  The  ultimate 

4 


test,  however,  is  comparison  with  in  situ  observations.  Equatorial  in 
situ  data  from  different  altitudes  will  also  shed  light  on  the  validity  of 
a  single  thin  screen  interpretation  of  an  integrated  measurement. 


Another  outstanding  issue  is  the  universal  applicability  of  the  two- 
component  spectral  model  refined  here  for  the  equator.  Since  the  physical 
processes  involved  in  determining  the  spectral  slopes  have  not  been  speci¬ 
fically  identified,  it  is  crucial  that  convective  structure  be  examined 


under  other  conditions;  e.g. ,  at  high  latitudes.  We  intend  to  elucidate 
these  issues  under  continuing  ONA  support. 
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